Tuning thermal transport in nanotubes with topological defects 
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Abstract 

Using the atomistic nonequilibrium Green's function, we find that thermal conductance of car- 
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strong Rayleigh scattering of high-frequency phonons. Phonon transmission across multiple defects 
behaves as a cascade scattering based with the random phase approximation. We elucidate that 
phonon scattering by structural defects is related to the spatial fluctuations of local vibrational 
density of states(LVDOS). An effective method of tuning thermal transport in low-dimensional 
systems through the modulation of LVDOS has been proposed. Our findings provide insights into 
experimentally controlling thermal transport in nanoscale devices. 



Carbon nanotubes (CNTs) are interesting materials for nanoscale electronic devices due 
to their outstanding electronic and thermal properties.-^ It is also found that CNTs have 
some unusual thermal properties, such as good thermal conductivity,-!^ thermal rectification- 
and thermal waveguide.™ Tuning thermal transport in low-dimensional systems like CNTs 
not only acts as an challenge for conventional nanoelectronics and energy conversion on 
a chip,-!^ but also gives rise to an emerging field such as phononic devices and thermal 
logic units.- Thermal conductivity has been tailored with some methods such as the disor- 
dered layer,-ii^ the isotopic mass disorder—"— and the pressure.— In contrast with the high- 
concentrations of isotopic mass disorder,—"— the structural lattice imperfection in CNTs 
can be well controlled^ during the growth or by the irradiation process. 

In this letter we aim to investigate the effects of multiple topological structure defects 
on thermal transport in CNTs through the approach of atomistic nonequilibrium Green's 
function(NEGF).— li^ We find that thermal conductance of CNTs can be efficiently tuned 
with only a few topological lattice defects. We propose an effective method of tuning thermal 
conductivity in low-dimensional systems through the modulation of local vibrational den- 
sity of states(LVDOS). Our findings are interesting to the experimental control of thermal 
transport in nano devices. 

We consider a carbon nanotube connected to heat reservoirs with different temperatures 
at each tube end. An energy current flows from high temperature to low temperature across 
the tube. To eliminate the interfacial thermal resistance between the tube and the reservoir, 
the heat baths are simulated with the same semi-inflnitely long tubes. Under steady state, 
the heat current / is described by NEGF.— Here we concentrate on the behavior of phonon 
scatterings induced by lattice imperfections so that nonlinear scatterings are ignored. In 
fact, nonlinear phonon interactions can be included through the calculation of the nonlinear 
self-energy.— If the elastic scattering is considered, the heat current expression is reduced 
to the Landauer formula with thermal conductance Gth at temperature T deflned as^ 

1 f^ df 

G.=-y^ d^fi^TM^. (1) 

where the phonon transmission T[uj] at a given frequency u can be calculated through the 
Caroli formula T = Tt^G^TlGTr). Here the self-energy density r^ = z(llj^ — T,'^), a = L, R 
for the left/right lead, C and C^ are the retarded/advanced Green's function, and f the 
Bose distribution for phonons. 



To calculate the NEGF atomistically, the dynamic matrix are derived from the second- 
order derivative of the second-generation Brenner potential^i^ with respect to displacements 
after the structures are optimized with the same potential. The Brenner potential can 
well reproduce the phonon dispersion relations for CNTs.^^ In comparison with the first- 
principle density-functional theory approach, this empirical potential is short-ranged and 
more tractable for systems with large number of atoms. ^^ During our simulations, there are 
thousands of atoms, which are beyond the ability of any first-principle method. 

The calculated phonon transmissions for the CNT (10, 10) through the NEGF are shown 
in Fig. [U^a) and Fig. Wljo)- With the increase of the number of structural defects, it can be 
seen from the figure that phonon transmission for high-frequency decrease rapidly, irrespec- 
tive of Stone- Wales or vacancy defects. To understand this behavior, we first investigate 
the phonon length of mean-free path(MFP). For perfect CNTs, the length of phonon MFP 
is estimated^ to be in the order of magnitude of micrometers at room temperature. Here 
the length of CNTs considered is a few nanometers. Thus, the transport process across the 
defect scattering can be described as ballistic-diffusive,-'^^'^^ where the total phonon trans- 
mission at a given frequency co can be written^'^'^ as T[u] = Tp[u]lo/{L + Iq). Here Tp is 
the phonon transmission for the perfect tube at frequency u, Iq the length of phonon MFP, 
and L the length of CNT. This formula has been used for estimating the length of phonon 
MFP.— The calculated length of phonon MFP is shown in Fig. [H^c), where we can find 
that the length of MFP decreases with the frequency in the quadratic form (oc w^^). For d 
dimensional systems, the Rayleigh scattering theory^^ gives the length of the MFP as a func- 
tion of oc co^^'^^^'. The CNTs can be assumed as quasi-one-dimensional systems with d = 1 
such that the phonon MFP has a quadratic dependence(cx w^^) by the theory of Rayleigh 
scattering. The frequency dependency of the phonon MFP in Fig. [T]^c) agrees well with 
the quadratic relation so that the phonon scattering by defects in CNTs can be understood 
in terms of Rayleigh scattering. We thus can conclude that the larger decrease of phonon 
transmission in the high-frequency range results from strong scatterings of high-frequency 
phonons, due to the rapid decay (oc w^^) of MFP characterized by Rayleigh scattering. 

Next we consider the behavior of phonon transmission across multiple topological lattice 
defects. In the ballistic-diffusive transport regime, the phonon transmission T/v across A^ 



defects can be described by the cascade scattering modeU^'i^ 

1 N N -I 
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where 7i is the transmission across single defect, N the number of defects, and Tp the 
transmission for the perfect tube. It is straightforward that the quantity {Tp — 7i)/7i in 
Eq. ([2]) has the additive property with the transmission probabihty 7i for each scatter 
placed in cascade. ^^ The calculated phonon transmission from the cascade model is shown 
in Fig. [T](a) and Fig. [U^b) denoted by the green dashed lines. We can find that phonon 
transmissions derived from Eq. ([2]) agree well with the atomistic calculations by NEGF for 
both Stone- Wales and vacancy defects. After performing computations for longer CNTs 
with more structural defects, we have also observed such consistency between the direct 
numerical calculations by NEGF and the predictions of the cascade scattering model. The 
reason for the good validity of the cascade model can be explained as follows. Note that 
Eq. ([2]) is associated with a random phase approximation for wave interferences.— Unlike 
the point isotope-disorder, structural defects introduce a large area of lattice distortion such 
that it is difficult to keep the phonon waves coherent. Therefore, the random phase serves 
as a good approximation for the scattered waves. When the concentration of structural 
defects is not too high, we think that the cascade scattering model can well describe phonon 
scatterings over multiple structural defects. 

Further, we elucidate the mechanisms of the scattering of phonon by structural defects 
through calculating the phonon density of states(DOS) using NEGF. The phonon DOS 
and LDOS can be expressed as the imaginary part of the Green's function as DOS[a;] = 
— 2a;ImTr(G'')/7r and LDOS[tj];; = —2ujIniTTx.y,z{G^')i.i/'^, respectively. The subscript / is 
the index of the diagonal term of the Green's function matrix for each atom. The trace 
operation Tr^^y.z for each atom in the LDOS is only carried out in three directions. It can 
be seen from Fig. [2]^a) that there is an appearance of some extended spectra of vibrational 
states for the defective tubes at a; > 1748cm~^. These extended spectra correspond to the 
localized vibrations of the defect lattices. Fig. Mjo) shows a profile of the LDOS. We can find 
that the perfect CNT has a uniform distribution of LDOS, while the LDOS for the defective 
tubes fiuctuates with the positions of atoms, especially for atoms located in defects. The 
different phonon modes are scattered by the fiuctuation of LDOS as shown in Fig. [2](b). 
An atomistic view of LDOS is further demonstrated in Fig. [2t^c). The LDOS provides a 



direct view of spatial fluctuations of vibrational states around the defective tubes. The local 
density of states for electrons can be measured by scanning tunneling microscopy, but the 
experimental measurement of the phonon LDOS has not been realized. We think that the 
measurement of local vibrational density of states is not only important to understanding 
thermal transport in nanoscale systems but also is vital to developing nano devices. 

Finally, thermal conductance as a function of temperature is shown in Fig. |3l The solid 
line on the top represents the ballistic upper bound for the thermal conductance of the 
perfect CNT. With the increase of number of defects, we can find that thermal conductance 
rapidly decreases and varies little with temperature at high temperature. The reason for 
this large reduction and small variations is that high-frequency phonons which contribute 
to thermal transport at high temperature is strongly scattered by the structural defects in 
terms of Rayleigh scattering. In contrast to thermal conductance of CNTs with isotope 
disorders,— "—li^ we find that even a few structural defects in CNTs can lead to a strong 
suppression of thermal transport by one order of magnitude. We further propose that the 
topological structural defects can offer an effective method of tuning thermal transport in 
low- dimensional systems like CNTs through the modulation of LDOS. One advantage of 
such approach is that topological lattice defects can be well controlled during the growth 
or by irradiations. Compared with the molecular dynamical simulations,— "2^ a decrease of 
thermal conductance with temperature is not observed at high temperatures due to the 
neglect of nonlinear phonon-phonon interaction during our calculations. 

In summary, using the atomistic NEGF we find that thermal conductance of CNTs 
with presence of topological imperfects is remarkably reduced. High-frequency phonons 
are strongly scattered by structural defects, which can be characterized by one- dimensional 
Rayleigh scattering. In comparison with that of single structural defect in CNTs,—""— we 
show that phonon scattering across multiple defects can be well described by the cascade 
model with the random phase approximation. We also demonstrate that phonon scattering 
is related to the spatial fluctuations of LDOS. An effective method of tuning thermal trans- 
port in low-dimensional systems by modulating LDOS has been proposed. We expect our 
findings can contribute insights to the experimental control of thermal transport in nanoscale 
devices. 
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FIG. 1: (Color online). Figures (a) and (b). Phonon transmission as a function of frequency for the 
defective CNT (10, 10). The green dashed lines in figures are derived from the cascade scattering 
model. Figure(c). The frequency dependence of the length of phonon mean free path(MFP) for the 
CNT with one Stone- Wales defect. The red dashed line illustrates the quadratic relation (oc OJ^'^) 
for a guide view. 
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FIG. 2: (Color online). Figure(a). The phonon density of states(DOS) as a function of frequency. 
Figure(b). A profile of local vibrational density of states(LVDOS) at frequency w = 236.9Gcm~^ . 
Figure(c). Atomistic view of LVDOS for the CNT(10, 10) with a SW defect and a vacancy defect 
at frequency uj = 1^23.12071"^ , respectively. 
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FIG. 3: (Color online). The temperature dependence of thermal conductance for the perfect and 
defective CNTs. The number of defects are indicated by the number in the circle. The color dashed 
lines represent thermal conductance for the CNTs with vacancies while the color solid lines are for 
the SW defective tubes. 
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